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Abstract — In this contribution, the 3-dimensional (3D) channel 
characteristics, particularly in the elevation domains, are ex- 
tracted through measurements in typical urban macro and micro 
environments in Xi'an China. Stochastic channel model param- 
eters are obtained based on the high-resolution multi-path pa- 
rameter estimates. In addition, a modified spatial channel model 
(SCM) for the line-of-sight (LoS) scenario is proposed where 
the LoS polarization matrix is parameterized in accordance 
with the reality. Measurement results justify the reasonability 
of the proposed model. These works significantly improve the 
applicability of the ITU SCM models in realistic 3D channel 
simulations. 
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Figure 1. The system diagram of the measurement set-up 
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Index Terms — 3-dimensional channel models, elevation of arrival, 
elevation of departure, line-of-sight, and polarization matrix 

I. Introduction 

. Multiple-input Multiple-ouput (MIMO) enhancements are pro- 
posed in the 3rd Generation Partnership Project (3GPP) Re- 
leases 10 and 11 to support eNodeB (eNB) antenna configu- 
■ rations capable of adaptation in azimuth. Recently researches 
in enhancing system performance through the use of antenna 
systems having a two-dimensional array structure has been 
paid lots of attention in order to provide more spatial degree of 
freedom in both elevation and azimuth domains 1 1 1. Additional 
control in the elevation dimension enables new transmission 
'techniques by using the 3-dimensional(3D)-beamforming or 
with the massive MIMO, such as the sector- specific elevation 
beamforming (e.g., adaptive control over the vertical pattern 
^ beam width and/or downtilt), advanced sectorization in the 
vertical domain, and user-specific elevation beamforming tech- 
niques. In order to evaluate the performance of the systems 
' utilizing these techniques, new channel models characterizing 
' channels in both vertical and horizontal dimensions in various 
environments with diff'erent user locations are necessary. 

Existing directional channel models as proposed in 3 GPP 
TR25.996 d, the IMT- Advanced standards and in the 
WINNER reports p] have been widely adopted for devel- 
opment of wireless communication systems. However, with 
respect to the elevation characteristics of the radio channel, 
these models are considered to be insufl&cient. For example, 
in the WIINNER-h project |5|, 3D channel modeling method- 
ology and parameters were proposed based on literature survey 
rather than extensive measurements. Thus, the applicability of 
the models suggested are questionable due to some obvious 
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errors, such as that the cross-correlation matrix of large scale 
parameters (LSPs) is non-positive definite. The elevation char- 
acteristics, such as mean elevation angular spread (MEAS), 
extracted based on measurements have also been reported 
in literature [T|, However, these models focus on 
the marginal elevation characteristics, rather than the joint 
characteristics in both azimuth and elevation, and thus, are 
lack of applicability for 3D MIMO system-level evaluation. 

In this contribution, 3D channel measurement campaigns con- 
ducted by Huawei Propagation Research Group are described, 
which have been performed in diff'erent urban scenarios with 
the objective of extracting the full 3D channel characteristics 
particularly at the eNB. New measurement results for the 
statistics of the LSPs, and the cross-correlation matrix thereof 
are elaborated. Furthermore, a questionable setting specified 
in the WINNER-h SCM models for the polarization matrix 
of the line-of-sight (LoS) path is corrected and verified with 
measurement data. 

The rest of the paper is organized as follows. Section |II] 
introduces the measurement setup and the environments where 
the measurements were carried out. In Section |IV] the channel 
characteristics extracted for the urban macro (UMA) and urban 
micro (UMI) scenarios are described and the updated values 
for the ITU IMT-A model parameters are reported. In Section 
|Vl a mistake in the SCM, WINNER-F, and ITU IMT- Advanced 
models for the LoS scenario is discussed, and the correction is 
proposed. Finally conclusive remarks are provided in Section 
lYll 



II. Measurement setup and environments 

A diagram for the measurement setup is illustrated in Figure 
[T] where the Agilent E4438C signal generator is applied 
as an eNB, generating a pseudo-noise signal with 35MHz 
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(a) Measurement routes in Area A (Picture: © Baidu Maps) 



Figure 2. Antenna arrays at eNB (left) and UE (right) sides 

bandwidth and center frequency of 2.6 GHz. Figure [2] depicts 
the photographs of the antennas used in the eNB and the user 
equipment (UE). The eNB was equipped with a planar antenna 
array with 32 antenna elements configured as 16 patches, each 
consisting of two antennas with ±45° polarizations, 7 dBi gain 
and 90° beam width on both horizontal and vertical planes. A 
crown-shaped antenna array is used in the UE, which has 50 
antenna elements allocated on 25 patches with the structure 
similar with those on the eNB antenna array. The horizontal 
and vertical spacings between the nearest neighboring antennas 
are 0.5 wavelength. The down- tilting angle of the eNB antenna 
array is approximately 7°. 

The measurement was conducted emulating the downlink 
scenarios, i.e. the eNB and the UE were considered as 
the transmitter (Tx) and the receiver (Rx) respectively. The 
high-resolution channel parameter estimation algorithm Space- 
Alternating Generalized Expectation-maximization (SAGE) 
algorithm ||9l was applied to estimation of 14 parameters for 
individual paths, i.e. the delay, Doppler frequency, direction 
(azimuth and elevation) of departure (DoD), direction of 
arrival (Do A) and polarization matrix. 

The measurements were conducted in the "High-tech" district 
of Xi'an, China. The UE was moving along two routes in the 
area A and area B as indicated in Figure [3] and |4] respectively. 
The photographs taken in the eNB side for the area A and 
B are shown in Figure [3] and |4] It can be observed that the 
area A and B are close to the definition of UMA and UMI 
scenarios specified in the WINNER modeling lH receptively. 
The heights of the eNB antenna in area A and area B were 
40 m and 14 m above the ground respectively. Notice that 
considering the limited eff'ective radiation of the antenna array 
in the eNB, the estimation range of azimuth of departure 
(AoD) and elevation of departure (EoD) are confined to be 
[-180°, 0°] and [0°, 180°], respectively. Furthermore, the AoD 
and EoD of the direction normal to the array plane are -90° 
and 90° respectively. 

III. Verification of the parameter estimation 

To ensure that the path parameters are correctly estimated, 
the angular estimates obtained in the LoS scenario in the 
case 4 in the area B are compared with their counterparts 
calculated geometrically based on the map information and 
the antenna heights. It can be observed from Figure |4] that 




Figure 3. Measurement routes and a photograph taken in Area A 
(a)^ Measurement routes in Area B (Picture: © Baidu Maps) 




Figure 4. Measurement routes and a photograph taken in Area B 
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(a) Estimated AoD of the LoS path 
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(b) Estimated EoD of the LoS path 
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Figure 5. Estimated AoD and EoD for LOS path in case 4 



in the measurement case 4, the UE started from the location 
A, moving to the location B, then back to A, and finally 
arrived at the location C. The theoretical AoD and EoD for the 
LoS path at the start and the end point of this measurement 
route are illustrated in Figure |4] Figure |5] (a) depicts the AoD 
estimates for the LoS path, which is selected to be the strongest 
path estimated by using the SAGE algorithm, versus the 
measurement snapshots. It can be seen that the AoD estimate 
increases gradually from -85° to -98°, and decreases to -85° 
when the UE arrived at the location A, and then increase 
again up to -70° when the UE arrived at the location C. This 
variation is consistent with the AoD's trajectory calculated 
theoretically based on the UE displacement. It can also be 
observed from Figure [5] (a) that there is a sudden change 
of the AoD at the 180th measurement snapshot. This is due 
to the fact that the orientation of the Tx antenna array was 
adjusted manually at that moment in order to make sure that 
the UE can be covered by the array's main lobe. The array was 
tuned back to its original position finally. Figure [5] (b) depicts 
the estimated EoDs of the LoS path versus the snapshots, 
which also fits well with the theoretically calculated EoDs. In 
addition, it is worth mentioning that the estimated trajectories 
of AoDs and EoDs are not smooth due to the limited intrinsic 
resolution of the measurement equipment in both delay and 
angular domains. 

IV. Measurement results 

A. Measurement-based 3D channel parameters 

The measurements were conducted with the objective of 
extending the ITU channel model 1 1 1 to include the elevation 
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(a) Urban macro LoS scenario 
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(b) Urban micro LoS scenario 
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Figure 6. EoD distribution in UMA & UMI LOS scenarios 

Table I 

Channel model parameters extracted 



Scenarios 




UMI 


UMA 






LoS 


NLoS 


LoS 


NLoS 


Delay spread 
(DS) logio([^]) 


n 

(T 


-7.08 
0.34 


-7.6 
0.8 


-6.66 
0.56 


-6.51 
0.2313 


ASD logio([°]) 


fi 
(T 


0.84 
0.19 


1.1 
0.3 


0.61 
0.32 


1.05 
0.1 


ASA logio([°]) 


fi 
(T 


1.64 
0.17 


1.5 
0.13 


1.68 
0.18 


1.87 
0.11 


ESD logio(n) 


fi 
(T 


0.27 
0.29 


0.6 
0.2 


0.11 
0.2 


0.8346 
0.0452 


ESA logio(n) 


fi 
(T 


1.06 
0.18 


0.88 
0.16 


1.49 
0.43 


1.26 
0.16 


Mean EoD ([°]) 


-2 


-2 


2 


-2 


Mean EoA ([°]) 


2 


2 


6 


10 


/T-factor [dB] 


fl 
(T 


9 
5 


N/A 
N/A 


9 

3.5 


N/A 
N/A 


EoA, EoD distribution 


Laplacian 


Cluster ESD 


1.45 


3 


1.16 


2.332 


Cluster ESA 


6.68 


7 


15.5 


7 


Per cluster shadowing std 


3 dB 


3 dB 


3 dB 


3 dB 



ASD: Azimuth Spread of Departure 
ASA: Azimuth Spread of Arrival 
ESD: Elevation Spread of Departure 
ESA: Elevation Spread of Arrival 
EoD: Elevation of Departure 
EoA: Elevation of Arrival 



dimension in order to generate the 3D channel model for 
system level simulations. Table [I reports the statistical channel 
parameters extracted for both the elevation and azimuth do- 
mains at the eNB and UE sides, as well as the delay domains. 
In the table, DS, ASD and ASA represent the delay spread, az- 
imuth spread for departure and for arrival, respectively. These 
parameters can be used to update the corresponding entries of 
the table A 1-7 in O for ITU channel model standards for the 
environments considered in the measurements. 
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Table II 

Mean azimuth angular spread (MA AS) and mean elevation angular spread 

(MEAS) AT THE eNB 





MA AS [°] 


MEAS n 


UMA LoS 


4.6 


1.77 


UMA NLoS 


11.2 


6.8 


UMI LoS 


11.54 


2.97 


UMI NLoS 


13.3 


6.28 



The 3D MIMO performance evaluation also requires the power 
elevation spectrum and elevation spread at the eNB. Table HI] 
reports the mean azimuth angular spread (MA AS) and mean 
elevation angular spread (MEAS) at the eNB. Figure [6] depicts 
the EoA distributions for the UMA and UMI LoS scenarios. It 
can be observed from Figure |6]that the Laplacian distributions 
can be used to fit the empirical distribution derived from 
measurements. 

Table [Till reports the cross-correlation coefficients between 
diff'erent large-scale parameters (LSPs) based on the measure- 
ment results. Here, the cross-correlation coefficient of two 
random variables, say a and b, are calculated as 

Y,(ai - a)(bi - by 

C(a,b) = 

where a and b denote the mean of a and b respectively. The 
3D channel realizations can be generated using the parameters 
given in Tables U to [Till following the ITU channel generation 
procedure. 

V. Generating channel coefficient for LOS path 

In the case where the elevation domain is taken into account, 
the contribution of the ^th cluster in the narrowband coeffi- 
cients of the channel between the ^-th Tx antenna and the uth 
Rx antenna can be calculated as lU 




+ <5(«-l)<t(0, (1) 



where H'^^^ Jj) represents the NLoS component and can be 
written as Q, the same as in |4, EQ.(4.20)] in the case where 
elevation domains are considered. Due to the limitation of 
space, the readers of interest are referred to [4 for detailed 
explanation of the notations arising in (0. 

In the case of clear LoS without obstructions, H^^J^^if) is 
written as ©, where the diagonal elements of the polarization 
matrix are co-phased, i.e. 

CoS = <oS='^LoS, 

which is diff'erent from the settings in the existing ITU and 
3 GPP geometry-based channel models that 0J^o^ and (p^^^ are 
generated independently. This revision is important since in 
the case where the LoS path is of free-space propagation, the 
polarization status of the transmitted EM wave shall be unaf- 
fected, which occurs if and only if the diagonal elements of 
the LoS path polarization matrix do not introduce extra phase 



Table III 

Cross-correlation of LSPs (DS,ASD,ASA, ESD and ESA ,SF,K) 



Scenarios 


UMI 


UMA 




LoS 


NLoS 


LoS 


NLoS 




ASD vs DS 


0.42 





0.21 


0.4 




ASA vs DS 


0.11 


0.4 


0.25 


0.6 




ASA vs SF 


-0.4 


-0.4 


-0.5 







ASD vs SF 


-0.5 





-0.5 


-0.6 




DS vs SF 


-0.4 


-0.7 


-0.4 


-0.4 




ASF) vs ASA 


0.11 








0.4 




ASD vs K 


-0.2 


N/A 





N/A 




ASA vs K 


-0.3 


N/A 


-0.2 


N/A 




DS vs K 


-0.7 


N/A 


-0.4 


N/A 


Cross- 
Correlations 


SF vs K 
ESD vs DS 
ESA vs DS 


0.5 

-0.07 

0.12 


N/A 
-0.5 





-0.3 
0.13 


N/A 
-0.3 





FSD vs ASD 





0.5 





0.3 




ESA vs ASD 





0.5 


0.44 


-0.2 




FSA vs ASA 


0.37 





0.3 







ESD vs ASA 


-0.26 





-0.2 







ESD vs SF 
















ESA vs SF 








-0.8 


-0.5 




ESD vs K 





N/A 





N/A 




ESA vs K 





N/A 





N/A 




ESD vs ESA 


-0.2 





-0.11 







DS 


7 


10 


30 


40 




ASD 


8 


10 


18 


50 


Correlation 


ASA 


8 


9 


15 


50 


distance 


ESD 


8 


9.5 


16.5 


50 


[m] 


ESA 


8 


9.5 


16.5 


50 




SF 


10 


13 


37 


50 




K 


15 


N/A 


12 


N/A 



SF: Shadow Fading 
DS: Delay Spread 



diff'erence to the vertical and horizontal polarized components 
being transmitted. The conventional setting of independent 
^LoS ^LoS ^'^sults in a consequence that the polarization 
status of received LoS component is diff'erent from that of 
the transmitted signal, which is inconsistent with the common 
experience in reality. 

To justify the amendment suggested, a measurement campaign 
was conducted to check the variability of polarization status 
for the signals received from the LoS path. Figure |7]depicts the 
premises of the measurement in the open rooftop environment 
with clear LoS between the Tx and Rx separated by 40 meters. 
The Tx was equipped with a patch consisting of two antennas 
transmitting ±45° polarized waves respectively. By feeding 
diff'erent phases into the input signals at radio frequency of the 
two Tx antennas, the transmitted wave may have desirable po- 
larizations, such as the linear, elliptical or circular polarizations 
with either left-hand or right-hand circulations. For example, 
the vertical linearly-polarized waves and circularly-polarized 
waves can be generated by tuning the phase diff'erence to 0° 
and 90°, respectively. Furthermore, during the measurements, a 
phasor-monitoring device was applied that makes sure the cor- 
rect phase diff'erence be maintained between the two branches 
of transmitted signals. The Rx was equipped with a dipole 
antenna, which has tunable orientation that allows receiving 
EM waves with any desirable linear polarizations. 

Figure [8] reports the received power at the output of the 
Rx with linear polarizations from 0° to 360°. Notice that 
the polarization slant angle of 0° is referred to the vertical 
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Figure 7. Measurement in a LoS environment. The antennas marked by the 
red and blue ellipses are for the Tx and Rx respectively. 

polarization direction in the measurement. It can be observed 
from Figure [8] that when the circularly-polarized wave was 
transmitted, the received power remained stable regardless of 
the polarizations of Rx antenna. When the vertical linearly - 
polarized wave was transmitted, i.e. the polarization was along 
0°, the received power achieved the maxima when the Rx 
polarization matched the Tx polarization, i.e. being 90° or 
270° polarized, and reduced to the minima when the Rx 
and Tx polarizations were orthogonal. It is evident that the 
polarization of the transmitted waves does not change in the 
LoS scenario considered in the measurement, implying that the 
diagonal elements of the LoS polarization matrix do introduce 
the same phase variation to both the vertical and horizontal 
polarized waves. 

VI. Conclusions 

In this contribution, the stochastic channel characteristics 
extracted based on measurements for urban macro and micro 
scenarios in Xi'an, China have been reported with emphasis in 
the elevation of arrival and of departure domains. The elevation 
angular spreads and the cross-correlation coefficients thereof 
with other large-scale channel parameters were calculated and 
compared with the existing 3 -dimensional (3D) ITU IMT- 
Advanced channel models for urban scenarios. The channels 
generated by using the updated model parameters yield non- 
negative definite correlation matrices which is difficult to 
obtain by using the existing models. In addition, modification 
were proposed for the conventional settings of the LoS- 
path polarization matrix defined in the WINNER-h/ITU/3GPP 
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Figure 8. Power of the received signal when the Rx polarization varies from 
0° to 360° polarizations 

spatial channel models. Measurement results justified the 
reasonability of the changes for real LoS cases. 
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